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1 INTRODUCTION 



In the present paper we will consider the relation between the deep inelastic 
structure function F2{x,Q'^) and the gluon density xG{x,Q^) in the region 
of small X. At very small x the main contribution to F2 comes from the 
sea quarks and the value of DIS cross section is given by the photon-gluon 
fusion subprocess Fig. 1. For large in the Leading Log Approximation 
(LLqA) and in a physical (say, axial) gauge it is enough to deal only with 
the diagram Fig. la. However within the LLog(-)A, which sums up all the 
terms of the type (a^ In ^)"', one has to take into account the graph Fig. lb 
also. 

Based on these graphs Fig. 1 one can extract the value of the gluon 
structure function G{x, k"^) from the experiment, where the DIS structure 
function ^2(0;, Q^) was measured. 

The question we are going to discuss here in more detail is: what are 
the values of gluon kinematical variables Xg and k"^ which correspond to the 
arguments x = xb and of the DIS cross section; i.e.F2? 

In Sect. 2 we will remind the LL^A formulae |1|] for the function F2 (and 
Fl) written in the framework of the high energy Ajj-factorization scheme [0], 
where the gluon k (see Fig. 1) is off mass-shell. 

Then by the numerical computation in Sect. 3 we will study the ratios 
of the deep inelastic (Bjorken) variables {xb and Q^) and the values of the 
gluon's momenta [xg and k"^ correspondingly), which are essential in the 
integrals of Eqs. (3) and (13); and how do these ratios depend on the form 
of the gluon structure function. 

To answer the last question we parametrize the gluon density by the power 
functions - xG{x, k"^) ~ x~^ ■ {k'^)"' and plot the logarithms of the ratios 

1 2 2 

(In 1^, In q2 , In qi) as the functions of powers A and 7. 

Even in the case of a strong scaling violation (7 = |) the average value 

< In p- > turns out to be rather large. It corresponds to k"^ ~ 0.8 2 GeV^ 
for = 16^64 GeV^. 

Thus there appears a danger that an important part of F2{x, Q^) is origi- 
nated from the region of a small gluon's transverse momentum (virtualities) 
{kt < 1 GeV), where the perturbative QCD approach is not justified. 

To follow the role of the infrared (small qt) region we have changed also 
the mass of a light quarkQThe value of F2 increases up to the factor of 1.5 -h 2 
when instead of the constituent mass rUq = 350 MeV one puts the small (close 
to the current) mass nig = 10 MeV. The situation becomes slightly better if 
we assume the saturation of gluon distributions at small x and fc^. 

^The heavy quark contribution to F2 was discussed in the previous paper [1]. 
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At very small x and fc^ the cross section a ~ ^'^^^2'^ ^ should tend to a 
constant. Thus at k"^ < (for example we put here /cq = 1 GeV^) the 
gluon density xG{x, k"^) ~ fc^; in other words at k"^ < IGeV"^ the anomalous 
dimension 7 = 1. In this case the rapid decrease of G{x,k'^) aX k"^ —>■ 
suppresses the role of the infrared region and the result becomes more stable 
(under the variation of m^). This point will be discussed also in Sect. 3. 

Finally we have to stress that the infrared {qt, kf < OA GeV) region gives 
not too small contribution to the ^Q^^^i^ ratio and one has to remember 
about it, trying to extract the gluon structure function from the F2{x,Q'^) 
experimental data. 



2 SMALL X DIS CROSS SECTIONS 

The cross section of heavy quarks as well as light quarks (if value is 
large enough) electroproduction is given schematically by the graphs in Fig. 
1. The main contribution to the cross section at small x is known to come 
from gluons. The lower ladder blocks presents the two-dimensional gluon 
distribution (p{x, k"^) which is a function of the fraction x of the longitudinal 
momentum of the initial hadron and the gluon virtuality. Its distribution 
over X and transverse momenta kx in hadron is given in the semihard theory 
[§] by function Lp{x,k'^). It differs from the usual gluon structure function 
G(x,g2): 

xG{x,q^) = / (p{x,ql)dql. (1) 



4^2 



TT 



Such definition of (f{x, k"^) makes possible to treat correctly the effects arising 
from gluons virtualities. The exact expression for this function can be ob- 
tained as a solution of the evolution equation which, contrary to the parton 
model case, is nonlinear due to interactions between the partons in small x 
region. 

In what follows we use Sudakov decomposition for the produced quarks' 
momenta pi^2 through the light-like momenta of colliding particles pa — Q-y 
and pb {j>a = Pb — ^) and transverse ones pi^2T'- 

Pl,2 = Xi^2PB + yi,2PA + Pl,2T ■ (2) 

The differential cross sections of a quark pair electroproduction^ have the 
^ o'ei,T,L denotes the virtual photon-proton transverse and longitudinal cross sections. 
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formQ: 

dael,T,L kernel ^ f ,2j .fj \ Xf , 1^ 

- i)2 J d kTd{kT - Pit - P2T) d{yi + y2 - I, 



dyldyl(PpiT(Pp2T (27r)4 (s) 

X ^4P^(A:^x)|Me,T,L^-. (3) 
fc2 y2 

Here s = 2paPb , k is the gluon transverse momentum and y* is the quarks 
rapidity in the 7*j9 c.m.s. frames, 

= e^S 1/2 = ^ rriT = v'pt + • 

|j^ei,T,Lp are the squares of the matrix elements |JlI| for the cases of transverse 
and longitudinal electroproduction, respectively. 
In LLA kinematic 

k ~ xpB + kx , (5) 

so 

,2 



k' ~ -4 . (6) 

(The more accurate relation is /c^ = —33^ but we are working in the kine- 
matics where x,y ^ 0). 

The matrix elements M are calculated in the Born order of QCD without 
standart simplifications of the parton model since in the small x domain 
there are no grounds for neglecting the transverse momenta of the gluon fc^ 
in comparison with the quark masses and transverse momenta, especially in 
the case of light quark pair production. In the axial gauge p%A^ = the 
gluon propagator takes the form D^^ik) = d^y{k)/k'^, 

d,M = ^^.v - Wb + ^""P^M^PBk) . (7) 

For the gluons in t-channel the main contribution comes from the so called 
'nonsense' polarization g^^^ which can be picked out by decomposing the 
numerator into longitudinal and transverse parts: 

= IWa + v'aPbMs + C ^ '^Ms ^ g;^. (8) 

The other contributions are suppressed by the powers of s. Since the sum 
of the diagrams in Fig. 1 are gauge invariant in the LLA, the transversality 



■^We put the argument of as to be equal to gluon virtuality, which is very close to the 
BLM scheme §; (see also g). 
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condition for the ends of gluon line enables one to replace by —k^rp/x in 
expression for g^,^. Thus we get 

d,^{k) ^ -2^ (9) 
X s 

or 

d,Ak) ^ 2^ (10) 
xys 

if we do such a trick for vector pb too. Both these equations for df^i, can 
be used but for the form @ one has to modify the gluon vertex slightly (to 
account for several ways of gluon emission - see Ref. ): 

Te// = — [i^ys - klr) Kt - klr^T + 2x (A;itA;2t)Pb]- (H) 
xys 

As a result the colliding gluons can be treated as aligned ones and their po- 
larization vectors are directed along the transverse momenta. Ultimately, the 
nontrivial azimuthal correlations must arise between the transverse momenta 
Pit and p2T of the produced quarks. 

In the case of photoproduction the alignment of the gluon makes the 
distribution of the difference of the transverse momenta of the quarks px = 
{pi — P2)t/2 to be proportional to 1 + cos'^9, where 9 is the angle between 
the vectors px and kx = (pi +P2)t/2. 

Although the situation considered here seems to be quite opposite to the 
parton model there is a certain limit in which our formulae can be trans- 
formed into parton model ones. Let us consider the case of pp collisions and 
assume now that the characteristic values of quarks' momenta pix and p2T 
are many times larger than the values of gluons' momenta, kix and k2T 

< Pit »< kiT> , < P2T »< k2T > , (12) 

and one can keep only lowest powers of kix, k2T- It means that we can put 
kiT = k2T = everywhere in the matrix element M except the vertices, as 
it was done in the parton model approach. 

On the other hand, the used assumption (|T2|) is not fulfilled in a more or 
less realistic case. The transverse momenta of produced light quarks should 
be of the order of gluon virtualities (QCD scale values). 
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3 RESULTS OF CALCULATIONS AND DIS- 
CUSSION 

The function /c^) is unknown at small values of k'^ . Therefore we rewrite 
the integrals over k2 in Eq. (|^) in the form 

r dkl5{kT - PIT - P2t) S{y, + ^2 - l)^p-^{e, x)\M\'- = 

Jo k^ 1/2 

AV27T'6ikT - Pit - P2t) 5(2/i + 2/2 - l)a.(gg)xG(x, Ql) f ^) - (13) 

V k-^ /fcT^o2/2 

+ r de^SikT - PIT - P2t) Kvi + 2/2 - l)^4?^(fc2, x)\M\'- , 
JQI k^ 2/2 

where Eq. ([1|) is used . 



Using Eq. (^) and Eq. (13) we perform the numerical calculation of light 
quark production. 

The results of computations are presented in Figs. 2-8. In Figs. 2 and 3 we 
plot the values of F2{x, Q^) coming from the sea u + u quark contribution in 
the cases of MT(S-DIS) P and GRV HO 0] gluon distributions, respectively. 
The parameter Qq in Eq. (13) which separates the region of 'small' and 'large' 
gluon momenta kt was chosen to be equal to 1 GeV (Figs. 2a and 3a) and 
2 GeV (Figs. 2b and 3b). It is seen that the results depend rather strongly 
on the value of Qq as well as on the mass of light quark and in both cases 
are differed from the original MT(S-DIS) (or GRV HO) predictions for u + u 
sea quark contribution (which, in our notations correspond to Qq > Q^ and 
rriq = 0) and shown by solid curves in Figs. 2 and 3. However the difference 
of our calculated F2{x,Q'^) values and the original ones is not large. It 
means that both the procedure of the gluon density extraction from F2{x, Q"^) 
experimental data and the accuracy of our calculations are reasonably enough 
and there is no a large contribution of any (say, nonperturbative) corrections. 

The only explanations is the contribution of the small transverse momenta 
kt < Qo and qt ~ ruq to the final value of F2. Even for Q^ = GAGeV^ the 
difference between the cases of rrig = 10 MeV and 100 MeV is of about 20% 
(up to 30% for GRV HO and Qq = IGeV) at x = 10"^ 

To understand better what are the typical quark (gj and gluon {kt,Xg) 
momenta, which are essential in the integrals Eq. (13), we have chose the 
parametrization of gluon distribution in the form 

X Y 



xG(a;,Q2) = i5 [ _1 — \ (14) 

^ ' VO.OOl/ \2QGeV^) ^ ' 
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and calculated the mean values of logarithms =< In ^ >, =< In ^ > 
and =< In ^ >. 

Xq 

The advantage of the logarithms is that they does not destroy the main 
power structure of the integrand (trying to estimate directly the value of 
< > or < p- >, one faces the divergence of the integral over dk'^ at large 
(small) k^ for the positive 7 < 1). 

The ratio of the sea w-quark to the gluon densities are given in Fig. 4 (in 
terms of ^2(0;, Q^) and xG{x,Q^) structure functions) for different values of 
7 (0.2 and 0.5) and A (0.2 and 0.5) in Eq. (14). At very small x < 3 ■ lO'^ 
this ratio tends to a constant and does not depend on x. However, even at 
large 7 = 1/2 (the value which corresponds to the asymptotics (extremum) 
of the BFKL solution) when from the formal point of view the convergence 
of the integral in the infrared region should be good enough, the dependence 
of the result on the light quark mass is still not negligible (~20% difference 
between the values of F2 for rrig = 350 MeV and 10 MeV at x ~ 10"'^ and 

= 16 GeV^). 

The situation becomes slightly better (see Fig. 5) if one assumes that for 
small the gluon density falls down as xG{x, Q"^) ~ (say, 7 = 1 for 

< 1 GeV^). Such a behaviour one might expect in the case of the gluon 
density saturation P] (i.e. a (Q') ~ ^-^^ ^ const at ^ 0). 

Unfortunately, the saturation hypothesis is also not enough to neglect 
completely the small qt contribution. Increasing the quark mass from 10 MeV 
to 350 MeV one diminishs the cross section more than 15% (for example, at 
X = 10-^ Q2 = 16 GeV2, A = 7 = 0.5; and 20% for 7 = 0.5 , without the 
saturation). 

This fact is confirmed by the 7-dependence of the function Vk (Fig. 7a). 
As it should be expected the value of rk depends slowly on A but noticiably 
increases with 7. For a more heavy photon (Q^ = 64 GeV^) becomes 
smaller. It means that the typical values of kf increases with not as 
quickly as (more close to y/Q^). In the best case (A = 7 = 0.5 for > 
1 GeV2 and 7 = 1 for < 1 GeV^) the value of (^ -1.7 at = 15 
GeV^) corresponds to k^, essential in the integral, of about 3 GeV^ (smaller 
than the usual cutoff Ql = 4-^5 GeV^ for the MT or MSR parametrizations 
[|, H of the structure functions). Without the saturation at = 16 GeV^ 
the rfc -3; i.e. typical k'^ ~ 0.8 GeV^. 

In order to estimate the values of the longitudinal gluon momenta (xg) we 
plot in Fig. 8 the logarithms of the xb to Xg ratios. The last ratios increase 
with as for a small the transverse momenta of quarks sometimes 
becomes larger than the photon virtuality {\/Q^), see Fig. 6, and one needs 
a large Xg of gluon to put these quarks on mass shell. In the region of ~ 
16 GeV^ the value ~ —1-2; it means that as an average the gluon Xg is 3 
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times larger than the Bjorken xb- This result slightly differs from the Prytz 
estimation where a gluon has Xg = 2xb- 

The value of F2{x,Q'^) is the sum of transverse {F2t) and longitudinal 
(-^2l) contributions. As one can see from the Table, at large the last 
contribution is more ifrared stable. 



4 CONCLUSION 

We hope that the ratio of F2 to xG{x, Q^) presented in Figs. 4, 5 as a 
functions of anomalous dimention 7 and intercept A may be used in future in 
order to extract, for example, the values of gluon density directly from the 
F2 structure function (in this case the values of 7 and A are also come from 
the F2 data as 7 = d In F2 /c? In and A = din F2/ dln{l/ x) . 

The results we have presented here are close to the standard gluon den- 
sities extracted from the derivative dF2{x,Q'^)/dln{Q'^) (using the Prytz 
method [Q, or the modified one, taken into account the NLO corrections 
[|10[] ). Nevertheless, here we better understand the typical values of gluon 
momenta we deal with and can control the role of the small kr (infrared) 
region. 

The mean logarithms Vg, and give us an estimate of the quark trans- 
verse momenta which are essential in the process, and the values of gluon 
virtuality (fc^) and proton momentum fraction [xg) which correspond to the 
point Xb and for the structure function F2. 

A lesson one has to extract from these calculations is the fact that an 
essential values of transverse momenta {qt) which do contribute to the next- 
to- leading order (NLO) corrections may be quite differ from the initial Q^. 
At not too large the value of may be even larger than (in Fig. 6 
rg > for Q2 < ig GeV^. 

However the main danger comes from the cases where the virtuality in the 
internal loops is much smaller than Q^. For example when (as it was done 
in [^) one includes in the anomalous dimension 7(0;) the NLO corrections 
of the type (as/u;)" in the first non-zero (4-loops) gluon-gluon correction ~ 
(as/w)^, he touches the region of the internal <C Q^; down to q^ < Q^/100 
due to the large numerical coefficient 5 = 14(^(3) in the Lipatov [12] diffusion- 
like distribution exp[— ln^(g^/Q^)/(45|/)] = exp[— ln^(g^/Q^)/(a • n)] where 
a = 6/ln2 ^ 24 and n is the number of loops. After this we use the value 
of the anomalous dimension 7(0;) (including the NLO corrections) as a local 
(from the point of view of InQ^ space) contribution at the given point Q^. 

There are no problems if is very large but one can not to do this for not 
too large where the essential transverse momenta of intermidiate gluons 
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may be rather small: close to the infrared region (Q/IO ~ ^^qcd)- 
This work is supported by the Grand INTAS-93-0079. 
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Table 1 



Transverse and longitudinal contributions to F2{x, Q^) for uu pair production 
at several points. 





X 


rrig, MeV 


= 4 GeV" 


= 64 GeV" 










F2L 


F2T 


F2L 


Without 


0.0016 


10 


0.126 


0.0354 


.392 


0.105 






350 


0.0816 


0.0242 


.346 


0.102 


saturation 


0.0064 


10 


0.0593 


0.0170 


0.186 


0.505 






350 


0.0378 


0.0116 


0.163 


0.489 


With 


0.0016 


10 


0.102 


0.331 


0.369 


0.104 






350 


0.0763 


0.0230 


0.338 


0.101 


saturation 


0.0064 


10 


0.0479 


0.0159 


0.175 


0.0500 






350 


0.0352 


0.0110 


0.159 


0.0484 
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Figure captions 



Fig. 1. Low order QCD diagrams for quark pair electroproduction via 
photon-gluon fusion. 

Fig. 2. u + u sea quark contributions to F2{x, Q'^) as given by MT S-DIS 
[^] sea quark distributions (solid curves) and our calculations using MT S-DIS 
gluon distributions for irtq = 10 MeV and 100 MeV (dashed and dash-dotted 
curves, respectively). The values of Qo parameter equal to 1 GeV (a) and 2 
GeV (b) were used in Eq. (13) 

Fig. 3. The same as Fig. 2 but for GRV HO parton distribution. 

Fig. 4. The ratios F2{x, Q^) {u + u sea quark contribution) to xG{x, Q"^) 
calculated from Eqs. (13), (14) for different values 7 and A. 

Fig. 5. The same as Fig. 4 assuming that at < 1 GeV^ 7 = 1 (i.e. 
saturation effect). 

2 

Fig. 6. Tq =< In^ > dependences on 7 (a) and A (b). 

Fig. 7. Tfc =< In^ > dependences on 7 (a) and A (b). 
Fig. 8. Ta, =< In— > dependences on 7 (a) and A (b). 
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Fig. 2b 



X 




Fig. 3a 



X 




Fig. 3b 



X 



F2/xG(x,Q^) for light quarks, 7*p, 



o 

X 



7 - 0.2 

A = 0.2 



a 

X 




-2 







3 
2 
1 



q2 = 64 GeV' 



1 6 GeV" 



Q' = 4 GeV' 



m,^ = 10 MeV 
m = 150 MeV 



r^, = 350 MeV 







10 



-5 



10 



-4 







-3 



10 



-2 



Fig. 4a 
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Fig. 4b 
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Fig. 4c 
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Fig. 4d 
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Fig. 5a 
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Fig. 5b 
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Fig. 6a 
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Fig. 8a 
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